One-particle interchain hopping in coupled Hubbard chains 
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Interchain hopping in systems of coupled chains of correlated electrons is investigated by exact 
diagonalizations and Quantum-Monte-Carlo methods. For two weakly coupled Hubbard chains 
at commensurate densities (e.g. n=l/3) the splitting at the Fermi level between bonding and 
antibonding bands is strongly reduced (but not suppressed) by repulsive interactions extending to 
a few lattice spacings. The magnitude of this reduction is directly connected to the exponent a of 
the ID Luttinger liquid. However, we show that the incoherent part of the single particle spectral 
function is much less affected by the interchain coupling. This suggests that incoherent interchain 
hopping could occur for intermediate a values. 
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The possible occurence of incoherent single-particle 
hopping in strongly correlated systems is presently ac- 
tively debated. Such an unusual behavior could, for 
example, explain the unconventional c-axis conductivity 
found in two-layered high-Tc compounds . In addition, 
the high transition temperatures of these materials may 
also be connected to incoherent hopping ||]: as a conse- 
quence of incoherence, the coupled planes can lower the 
total energy by coherent pair-hopping processes, leading 
to an amplification of the tendency towards supercon- 
ducting pairing. 

Perturbation and renormalisation-group treatments 
[^-^ of the transverse hopping between chains reveals 
that the key parameter which governs interchain hop- 
ping is the non-universal exponent a of the ID Luttinger 
liquid (LL) |^ characterizing the low frequency behav- 
ior of the density of state N{u)) ~ w". For a > 1, ie 
for sufficiently strong interaction, the interchain hopping 
becomes irrelevant. Recently Clarke, Strong and Ander- 
son suggested a somewhat different criterium: let us 
consider, for (real) times r < 0, a system of two de- 
coupled chains with different electron numbers on each 
chain in some initial state |5'o)- Then, if the interchain 
hopping t±^ is switched on at r = two different behav- 
iors can be observed: the quantity P{t) = |yl(T)p, with 
■^i''') = (^o|e'^^|4'o) corresponding to the probability 
for the system to be in its initial state, at a later time 
T > might show oscillatory or monotoneous behav- 
iors characterizing "coherent" or "incoherent" interchain 
hopping, respectively. It was suggested that such an in- 
coherent behavior could occur at values of a as small as 
0.5. 

The first numerical attempt to examine the question 
of coherent versus incoherent single particle hopping was 
initiated by two of us By exact diagonalizations 

of two coupled t-J chains, an interesting connection be- 
tween coherence and integrability was found. However, it 
should be emphasized that, in this approach, the initial 
state l^o) at T = is far from equilibrium. Indeed, con- 



structing the whole system with t± = Q corresponds to 
a macroscopic perturbation for the full hamiltonian (in- 
cluding the tranverse hopping). Interchain single particle 
hopping in the true GS (tj_ ^ 0) is still an open question. 

For this purpose, let us first consider two weakly cou- 
pled chains described by an extended Hubbard model, 



Vrtli^Xni+r.X ~ *i 5Z (4,l,aCi^2,o- + h.c), (1) 



where c\^^\s& creation fermion operator at site i of the 
chain A (=1,2). We shall deal here with two cases; (i) 
a purely on-site interaction (V^ — 0) or (ii) an extended 
repulsion of the form Vr = U /{r + 1) limited to r=l and 
2. 

In contrast to previous numerical work, we shall as- 
sume that, at r < 0, the system contains M electrons 
in the ground state (GS) |0o) of the full hamiltonian (|l|) 
and that, at r = 0, only a single extra particle (e.g. 
a hole) of longitudinal momentum k is added to, let 
say, chain number 1. Provided the occupation number 
nks = (0o|Cfe,i,cr'^fc,i:<Tl'^o) is non zero, the correspond- 
ing probability P{t) is then defined by the amplitude 
= /^oo^(^)e'""^^/'^fe.i> with 
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A{lo) = /m{(0o 



-fc,i. 



ie-H 



\^o)}- (2) 



This corresponds to the particular choice of an initial 
(normalized) state of the form |4'o) = ^^y/2 <^k.i,a\4'o}- 
Note that when an electron is added instead of a hole 
the creation and annihilation operators in (|^) have to be 
interchanged. At finite temperature T — 1/ f3 all expecta- 
tion values in |0o) must be replaced by thermal averages. 

Writing Ck_i^a in terms of the bonding and anti-bonding 
operators Ck.k±,rT {k± = 0, tt) it is straightforward to 
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express (g) as a linear combination, A{uj) = {A{k, 0, lj) + 
y4(/c,7r,a;))/(l + e^("~^)) where y4(k,w) = -^ImG{k,uj) 
corresponds to the usual spectral function in momentum 
space (k = {k,k±)). At this point, it is instructive to 
consider the non- interacting U=0 limit (and T — 0). In 
this case, A(k, oj) ~ d{uj — ek) and the two bands ek for 
= and TT are split by 2t±. If the momentum k 
(close to kp) is chosen in such a way that both fcj^ = 
and k± = TT states in |0o) are occupied, P{t) acquires 
an oscillatory behavior with a characteristic time scale 

For a generic LL the quasiparticle (5-peak is in fact re- 
placed by a singularity |^ corresponding to a collective 
excitation (in fact split holon and spinon peaks). How- 
ever, it is clear that an oscillatory behavior of P(t) is only 
possible if t± splits the dispersion relation of the bond- 
ing and anti-bonding modes by a finite amount although 
this alone might not be sufficient. One possible scenario 
would be that the whole spectra for = or tt (i.e. 
the poles of the one particle Green functions) would be 
displaced by ±t± respectively. In this case, P(r) would 
exactly take the form cos^(tj_T)Pi_D(r) where Pid{t) is 
independent of t± and is directly related to the ID spec- 
tral function AiD{k,u!) of the LL. However, we expect 
that, in general, (i) the incoherent part of AiD{k,uj) 
(i.e the part not related to the collective mode) is less 
sensitive to the transverse hopping at energies lo ^ t± 
and (ii) the relative weight of this part increases with 
the strength and the range of the interaction within the 
chains. Since analytic approaches of such a problem are 
extremely difficult, we shall rely in the following on exact 
diagonalizations of finite clusters (where dynamical and 
time dependent quantities can be calculated exactly) and 
on Quantum Monte Carlo (QMC) data. 

As a preliminary study, the various LL parameters of 
the ID single chain (e.g. a, Kp, spin and charge veloc- 
ities, etc..) are determined. The parameters of the ID 
on-site Hubbard model are known from numerical cal- 
culations based on the Bethe Ansatz |l^]. In the case 
of the extended model these parameters are obtained by 
standard finite-size scaling of ED results. For more tech- 
nical details the reader can refer to Ref. |11|. Hereafter, 
we shall restrict ourselves to a fermion density of n=l/3 
where relatively large values of the exponent a (up to 
^ 1.78) can in principal be obtained JT2| . Our numeri- 
cal data show that the maximum value can be realized 
in the case of the extended model for U ^ 10.5. Above 
this value a gap appears in the single particle spectral 
density possibly signaling a metal-insulator transition of 
no interest for us in the present work. 

The result for the spectral function of the double Hub- 
bard chain deduced by analytic continuation of the QMC 
data by maximum entropy methods is shown in Fig. (|l|) 
for U/t = 8. Bonding and anti-bonding bands dispersing 
through the fermi level {lo — fi — 0) can be clearly seen. 
Although error bars are rather large |^3| a splitting AE 
of order 2t± can be measured. It is a priori puzzling that 
such a large U does not lead to a substantial reduction of 



the splitting. However, one has to remember that a re- 
mains very small for an on-site interaction, e.g. a ~ 0.07 
for U/t = 8. 

As a comparison, the dispersions of the two low-energy 
branches corresponding to sharp peaks in A{h, uj) for 
k± = and k± ~ t: and obtained by ED techniques 
are shown in Fig. (|^) for the same value of U, in qual- 
itative agreement with the QMC results. Note that the 
boundary conditions along the chain direction (periodic 
or anti-periodic) are always chosen in such a way that 
the GS |(/>o) has a closed shell configuration. It should be 
emphasized that the band splitting picture is only valid 
sufficiently close to the Fermi level /i and that, on finite 
clusters, different values are obtained immediately above 
(i.e. |fc| > kp) or below (i.e |fc| < kp) the chemical po- 
tentail fi. 

A careful study of the ED data reveals that AE', for 
sufficiently small t± , behaves like at± + bt^ , where a and 
b depend on the parameters and system size. The ratio in 
the limit of vanishing t± is then easily extracted and the 
results are displayed in Fig. (||a) for the Hubbard and 
the extended Hubbard models and various system size. 
Assuming finite size corrections of the form l/L^ [ p^ , 
where is the length of the chains, and considering the 
average between the splittings above and below ^, one 
can also estimate AE/2t± in the thermodynamic limit 
as seen in Fig. (^d). Clearly repulsion at intermediate 
distances is much more effective to reduce the splitting. 
However, the plot of the same data vs a instead of U in 
(^d) suggests that a alone seems to be the key parameter 
controling the splitting. Indeed, an extended interaction 
leads to larger a values and hence to smaller AE/2t±. 

The previous study suggests that, at least for a < 1, a 
tranverse hopping t± leads to a splitting of the bonding 
and anti-bonding branches in agreement with perturba- 
tive treatments of t± (|-||]. However, the expression for 
P(t) involves the whole spectral function A(k, oj) at all 
frequencies so that it is not clear whether a spliting of 
the low energy singularity alone will necessarily lead to a 
coherent interchain behavior. The spectral functions at 
= and k± = TT oi the on-site Hubbard model and for 
a momentum k close to fc^- are shown in Fig. (^a). Most 
of the spectral weight is located at the low energy sin- 
gularity which is clearly split by the transverse hopping 
t±. On the contrary, in the case of the extended Hub- 
bard model (with parameters such that a ~ 0.53) shown 
in Fig. (IJb), more spectral weight appears far from the 
Fermi level and is less affected by t±. 

The probability P(t) has been computed from the 
knowledge of the spectral functions shown in Figs. (^,b) 
and results arc shown in Fig. (^). For the two types of 
models wc observe both fast and slow oscillations which 
have different physical origins. The fast oscillations are 
due to the fine structure of the ID spectral function 
AiD(k, w) (i.e. for t± = 0); for example, the sharp 
peaks of Fig. (^) might be attributed to shadow bands 
p5[ . We expect that, in the thermodynamic limit, these 
peaks will broaden and hence the fast oscillations will 
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get strongly damped if not suppressed. More interest- 
ingly, the slow oscillations have a pseudo-periodicity of 
order it / /S.E and are due to the effect of t^. They re- 
veal some degree of interchain coherence. However, in 
the case of the extended Hubbard model with a ~ 0.53 
a larger damping can be seen. It is not clear yet whether 
these slow oscillations will completely desappear in the 
thermodynamic limit. 

The previous results suggest that a splitting of the dis- 
persion of the single-particle low-energy mode of the LL 
by the tranverse coupling alone does not necessarily lead 
to an oscillatory behavior of P{t) which, according to 
Ref. [0 , characterizes a coherent interchain hopping. In- 
deed, for a sufficiently repulsive intermediate range (e.g. 
up to distances of 2 or 3 lattice spacings) interaction the 
large incoherent part of the single-particle ID spectral 
function might play a very important role in preventing 
coherent interchain hopping. 

DP acknowledges support from the EEC Human 
Capital and Mobility program under Grant No. 
CHRX-CT930332. WH acknowledges EC Grant No. 
ERBCHRX-CT940438 and the Bavarian FORSUPRA 
Progam. We also thank IDRIS, Orsay (France), for allo- 
cation of GPU time on the G94 and C98 GRAY super- 
computers, and the HLRZ Jiilich and LRZ Miinchen for 
access to the CRAY Y-MP supercomputers. 



[1] S.L. Cooper et al., Phys. Rev. Lett. 70, 1533 (1993); Y. 

Maeno et al., Nature 372, 532 (1994). 
[2] P. W. Anderson, Phys. Rev. Lett. 67, 3844 (1991). 
[3] H. J. Schuiz, Int. J. Mod. Phys. B 5, 57 (1991). 
[4] M. Fabrizio, A. Parola and E. Tosatti, Phys. Rev. B 46 

3159 (1992); C. Castellani, C. di Castro and W. Metzner, 

Phys. Rev. Lett. 69, 1703 (1992); M. Fabrizio, Phys. Rev. 

B 48, 15838 (1993). 
[5] C. Bourbonnais and L. G. Caron, Int. J. Mod. Phys. B 

5, 11033 (1991); D. Boies, C. Bourbonnais, and A.-M.S. 

Tremblay, Phys. Rev. Lett. 74, 968 (1995). 
[6] F. D. M. Haldane, J. Phys. C 14, 2585 (1981); see also 

J. Solyom, Adv. Phys. 28, 209 (1979). 
[7] D. G. Clarke, S. P. Strong and P. W. Anderson, Phys. 

Rev. Lett. 72, 3218 (1994); see also S. Chakravarty, A. 

Sudb0, P. W. Anderson and S. P. Strong, Science 261, 

337 (1993); D. G. Clarke and S. P. Strong, Sissa/cond- 

mat9510168 preprint (1995). 
[8] F. Mila and D. Poilblanc, Phys. Rev. Lett. 76, 287 

(1996). 

[9] V. Meden and K. Schonhammer, Phys. Rev. B 46, 15753 
(1992); J. Voit, Phys. Rev. B 47, 6740 (1993); K. Penc, 
F. Mila and H. Shiba, Sissa/cond-mat9501085 preprint 
(1995). 

[10] H. J. Schuiz, Phys. Rev. Lett. 64 2831, (1990). 

[11] F. Mila and X. Zotos, Europhys. Lett. 24, 133 (1993). 

[12] H. J. Schuiz, "Correlated Electron Systems", p. 199, ed. 



V. J. Emery (World Scientific, Singapore, 1993); H. J. 
Schuiz, "Strongly Correlated Electronic Materials: The 
Los Alamos Symposium 1993", p. 187, ed. K. S. Bedell, 
Z. Wang, D. E. Meltzer, A. V. Balatsky, E. Abrahams 
(Addison- Wesley, Reading, Massachusetts, 1994). 

[13] More accurate determination of the splitting from QMC 
can be obtained from the relative shift of the occupation 
numbers Wk for k± — and k± = tt. See H. Endres et al. 
(in preparation). 

[14] 1/Lx finite size corrections are in fact clearly established 
in the case of coupled spinless fermion chains; See S. Cap- 
poni et al. (in preparation). 

[15] K. Penc, K. Hallberg, F. Mila and H. Shiba, Sissa/cond- 
mat9512152 preprint (1995). 

FIG. 1. QMC result for A(k, cj) vs k = {k, fcx) on a 2 x 24 
Hubbard model with U/t = 8, t±/t = 0.15 at n = 1/3 and an 
inverse temperature of /3t — 10. Part (a) and (b) are density 
plots for the k± = and k± = tt branches respectively with 
the darker shading corresponding to higher spectral weight, 
and the points with error bars showing the position of the 
peaks. 

FIG. 2. Exact diagonalization result for the lowest exci- 
tation energies vs k = (fc, k±) obtained on a 2 x 12 Hubbard 
model with U/t = 8, txA = 0.15 at n = 1/3. The fuU symbols 
correspond to the main bonding and antibonding branches. 
Other low energy poles in Ak,u; (with smaller spectral weight) 
are also indicated by open dots. The thin dotted lines corre- 
spond to the U = case and the dashed line to the chemical 
potential. 

FIG. 3. (a) Normalized splitting AE/2t± extrapolated to 
tx — » (see text) vs U obtained by exact diagonalizations 
of 2 X 6 (open symbols) and 2 x 12 (full symbols) systems 
with on-site (O-S) Hubbard or intermediate range (I-R) re- 
pulsive interactions (see text). Both momenta k > kp (-I-) 
and k < kp (-) in the vicinity oi kp have been considered as 
indicated on the plot; (b) same quantity plotted as a function 
of a. The points with error bars correspond to crude esti- 
mates for the thermodynamic limit N ~* oo (see text). Error 
bars are estimated from separate extrapolations of the (-I-) 
and (-) data. 



FIG. 4. Spectral function A(k, k±,ui) vs on a 2 x 12 
system with n = 1/3 filling and for a longitudinal momentum 
k = 7r/12 < kp- The full and dashed lines correspond to 
k± = n and k± — 0, respectively. Only the uj < n region 
is shown (occupied states), (a) and (b) correspond to an 
on-site U — S Hubbard repulsion and to an intermediate range 
repulsion, U = 6, Vi — 3 and V2 = 2, respectively. 
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FIG. 5. Probability P(r) for a 2 x 12 system at n = 1/3 
filling and an extra hole with a longitudinal momentum 
k — 7r/12 < kp- The dashed and full lines correspond to 
an on-site Hubbard repulsion U = 8 and to an intermediate 
range repulsion, with U = 6, Vi = 3 and V2 = 2, respectively. 
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(3t = 10, (n) = 1/3 
A^/2t2; = 0.943 ± 0.068 



-0.5 



2x12 





At - A 2x6 OS (+) 
V 2x6 OS (-) 
>^ X 2x6 IR (-) 

H h 2x6 IR (+) 

^ A 2x12 0S(-) 
T 2x12 0S(+) 
■ ■ 2x12 IR(-) 
2x12 IR(+) 



0.0 





O O 2x6 (+,-) 
♦ ♦ 2x12(+,-) 
■ ■ 2xinf 



0.0 
0.00 



0.25 



T 



o 



♦ 



(b) 



a 



0.50 



0.75 



2 




1.0 




